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Statistical hadronization is presented as mechanism for (strange) parti- 
cle production from a deconfined quark-gluon plasma (QGP) fireball. We 
first consider hadronic resonance production at RHIC as a test of the model. 
We present in detail how the hadrochemistry determines particle multi- 
plicities and in case of sudden hadronization allows investigation of QGP 
properties. A comparative study of strange hadron production at SPS and 
RHIC is presented. The energy dependence of physical observables shows 
regularities and a potential discontinuity in the low RHIC range, when 
comparing these different energy domains. Considering the energy scan 
program at CERN-SPS we show that the J-k^ discontinuity is a baryon 
density effect. 

PACS numbers: 12.38.Mh,24.10.Pa,25.75.-q 



1. Introduction 

1.1. Setting the stage 

It is believed today that a new state of matter has been formed in 
relativistic nuclear collisions at BNL-RHIC and at CERN-SPS, and perhaps 
even in 1800 GeV pp elementary interactions. The question of considerable 
interest is if this is the hot quark matter (quark-gluon plasma, QGP) state of 
matter. The paradigm of QGP originates in the quantum many body theory 
of quark matter ^ 12] , which lead on to the formal recognition within the 
framework of asymptotically free quantum-chromodynamics (QCD) that at 
a very high temperature perturbative quark matter state must exist 
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Further reasoning based on a study of the 'boiling state' of dense hadron 
gas within the scheme of Hagedorn's statistical bootstrap of hadronic matter 
has lead from a different direction to the consideration of the transition to 
a hadron substructure phase The present day lattice-QCD numerical 

simulations allows to evaluate more rigorously the expected properties of 
a equilibrium state of QGP, as we today call hot quark matter [HI El IH1 
in]. QGP is the equilibrium state of matter at high temperature and/or 
density. The question is if, in the short time such conditions are available 
in laboratory experiments, this new state of matter can be created. There 
is no valid first principles answer available today and we must address this 
issue experimentally. 

It is obviously very hard to probe experimentally the QGP phase which 
exists a very short time. It is not widely accepted how a reaction involving 
formation of the deconfined state can be distinguished decisively from one 
involving reactions between individual confined hadrons. We observe in 
the final state of a high energy interaction reaction, irrespective if, or not, 
deconfinement has been achieved, a very large number of hadronic particles. 

The phenomenon of J and jet quenching probes the early stage of the 
deconfined phase. Its relation to deconfinement has been argued away in 
terms of effects of matter density and changes in primary production re- 
action mechanisms. Plasma photon radiation is buried in the background 
of final state photons from, e.g., i:^ — > 77. Dileptons (that is in fact vir- 
tual photons) suffer from large experimental uncertainties, and a final state 
background which depends of new physics as well (hadron properties at 
finite temperature and density). 

Perhaps the most promising path to a convincing evidence for the forma- 
tion of the deconfined state in relativistic heavy ion collisions is the precision 
study of the production of soft hadrons. Insights into the nature of the hot 
hadronic matter fireball are gained from interpretation of an array of par- 
ticle yields and spectra. A precise model description allows to evaluate the 
global properties such as energy, entropy, strangeness content of all parti- 
cles produced. This opens up to further investigation the properties of the 
matter state that has hadronized. 

Supporting this approach is the study of several reaction systems, e.g., 
comparison of AA with NN and/or NA reactions. In this case, the change 
of particle yield when this contradicts the hadron rescattering model expec- 
tations is evidence of new physics. Normally one encounters in an A A 
experiment enhancement of expected yield comparing to A^A^ reactions. 
More recently, some have argued that the AA yield of (strange) hadrons 
is 'normal' and that the effect is due to a 'suppressed' A^A^ base yield. This 
argument does not change the fact that in any case the AA yield does not 
arise from a sequence of A^A^ reactions. 
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1.2. Survey of this work 

The model we employ arises from the observation that hadron produc- 
tion is well described in a very large range of yields solely by evaluating the 
phase space size. Originally proposed by Fermi in order to describe multi- 
particle production in high energy cosmic ray interactions lllj. in the 
past 40 year the Fermi model (i.e., statistical hadronization) has been tried 
successfully in many particle production environments, beginning with p-p 
collisions at various energies, and branching out to heavy ion as well as the 
elementary non-hadronic e"'"e~ reactions. 

We begin, in section |2l with a discussion of the role hadronic resonances 
play in testing the principles of statistical hadronization. We survey, in sub- 
section 12.11 the latest experimental hadronic resonance production results 
from RHIC. We next describe how particle yields and ratios can be obtained. 
We introduce, in subsection 12.21 the ideas of hadrochemistry and show, in 
subsection 12. HI how particle ratios determine hadronization parameters. We 
turn to the study of phase space densities in subsection 12.41 We show that 
the entropy rich QGP phase hadronizes into chemical non-equilibrium pion 
gas that can absorb the entropy excess. For a more detailed discussion we 
refer to our earlier reviews |121 1131 [Hj . 

We introduce strangeness as signature of deconfined QGP phase in sec- 
tion El We discuss the interest in strangeness both from experimental and 
theoretical perspective in subsection 13.11 We introduce, in subsection 13.21 
the two step strange hadron production process, in which the hadroniza- 
tion benefits from ample supply of strangeness already produced. This is 
the origin of multi strange hadron enhancement. We show, in subsection 
13.31 that the observed enhancement effect is genuine. In subsection 13.41 we 
discuss how the outward collective flow of matter influences the resulting 
fireball hadronization conditions and assists the sudden hadronization. 

Section ^ presents our findings about the hadronization at SPS and at 
RHIC. We define first the method of analysis in subsection 14.11 and follow 
this through with a discussion of RHIC fits in subsection 14.21 We then 
address the energy dependence of the statistical parameters characterizing 
the fireball source in subsection 14.31 An interesting result that follows is 
about the behavior of energy per baryon stopping. Our primary interest is 
the strangeness yield which we address in subsection 14.41 

Our study leads directly to the final question, is there a energy thresh- 
old for the onset of quark-gluon plasma formation, which we address in 
section [SJ On important element of our study is to identify observables 
which are independent of the different baryon densities reached. We look 
from this perspective at the kaon to pion ratio discontinuity reported at 
low end of SPS energy range. We show that this is a result of baryon den- 
sity effect in subsection 15.11 In subsection 15.21 we discuss a mechanism of 



4 



zakop03c printed on February 1, 2008 



kaon to pion ratio enhancement arising from fast transverse expansion, and 
present the maximum enhancement expected at LHC. A brief evaluation of 
our findings and their importance follows in subsection 15.31 



2. Statistical hadronization 

2.1. Hadron resonances 

The initial test of statistical hadronization approach to particle produc- 
tion is that within a particle 'family', particle yields with same valance quark 
content are in relation to each other thermally equilibrated. Thus the rela- 
tive yield of, e.g., K*{sq) and K{sq) is controlled only by the particle masses 
mj, statistical weights (degeneracy) gi and the hadronization temperature 
T. In the Boltzmann limit one has (star denotes the resonance): 

N* _ g*m*^K2{m*/T) 

W ~ gm?K2{m/T) ' ^ ' 

Validity of this relation implies insensitivity of the quantum matrix element 
governing the coalescence-fragmentation production of particles to intrinsic 
structure (parity, spin, isospin), and particle mass. The measurement of 
the relative yield of hadron resonances is a sensitive test of the statistical 
hadronization hypothesis. 

The experimental ratio |151 1161 [T7j : 

(K* + W) (K* + W) 

^ ' = 0.26 ± 0.03 ± 0.07, ^ ' = 0.20 ± 0.01 ± 0.03, 

2iv zK 

has been measured, left for most central collisions with ^JsNN = 130GeV, 
and right for most central 200 GeV reactions, in RHIC-130, and RHIC-200 
runs respectively (65-1-65 GeV and 100-1-100 GeV per nucleon head on col- 
lisions of two nuclear beams). The RHIC result naturally arises from the 
ratio Eq. allowing for resonance decays within the statistical hadroniza- 
tion model for the temperature T ~ 145 MeV. This is the favored solution 
of a fit which allows for chemical non-equilibrium ^Sj, see subsection 14.21 
However, the yield: 

A(1520) 

-^ — ^ = 0.022 ±0.01, 
A 

recorded in the RHIC-130 and RHIC-200 runs ^H^Uj, is two times smaller 
compared to statistical hadronization expectation. Statistical hadronization 
at higher temperature than T = 145 MeV does not resolve this discrepancy, 
but adds to it a K* /K" suppression. 
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These two ratios, seen together, are, at first, difficult to understand. 
Let us first note that an apparent production suppression may actually be 
detectability suppression. Namely, should the decay products of resonances 
rescatter on other particles after formation, their energies and momenta will 
change. Hence, not all produced resonances can be, in general, reconstructed 
from the decay products energies Ei and momenta pi by testing for the 
invariant mass: 

m*^ = {Ei + E2f -{pi+P2?. (2) 

The rescattering effect depletes more strongly the yields of shorter lived 
states. These decay sooner and are thus more within the dense matter 
system. If resonance yield is as predicted by statistical hadronization, the 
decay life span is expected to be larger then the hadronization time, e.g., 

he 

eTK* = — = 3.9 fm > crhad- 
i K* 

This assures that most K* decay after thermal freeze-out of decay prod- 
ucts [K, vr) in free space. On the other hand, we recall that the lifespan 
of A(1520) is about 2.3 times longer than that of K* . This means that 
contrary to the expectation derived from rescattering effect, it is the longer 
lived A(1520) state for which we find the signature to be significantly de- 
pleted as compared to what could be maximally seen given the statistical 
hadronization yield. 

A natural explanation of the A(1520) signature suppression is in matter 
modification of the free space F = 15.6ibl MeV decay width [25 . This width 
corresponds to cr = 12.6 fm. Thus for sudden hadronization, implied by the 
observed yield of K* ., A(1520) decay should occur well outside the reaction 
domain. However, if in medium width is at the level of 200 MeV, within a 
hadronization lifespan of 0.7fm/c half of the produced A(1520) yield would 
become unobservable due to rescattering of decay products. Such a reso- 
nance width of 200 MeV is natural for a 400 MeV excited hadronic resonance 
state, once presence of the medium removes superselection constraints. 

Another proposed explanation of the X*, A(1520) riddle is that there is 
K* yield enhancement by regeneration ^22]. This vr + -fC — > K* production 
mechanism is generally associated with a long lived hadron phase. Factor 
5 enhancement over statistical hadronization yield is needed so that given 
the life span difference, the apparent reduction of A(1520) by factor 2 is 
consistent with the observed high K* yield. We are assuming here that the 
regeneration of A(1520) by A^, K collisions is not significant. However, a 
large K* yield enhancement contradicts basic principles, in that the heavy 
K* would compete with K in the overall bare yield. For this reason, we do 
not believe that this is a viable explanation. 
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We conclude that the statistical hadronization is doing reasonably well, 
the predictions in general agree with experiment in considerable detail, and 
where strong deviation is observed (A(1520)) a good reason for this is at 
hand. Resonances test both statistical hadronization, and the nature of the 
hadronization process |23]) and open a path to study in medium modifica- 
tions of particle hadronic widths. 

2.2. Chemical fug acities 

The freeze out temperature determines the shape of initially produced 
particle spectra, and thus it has considerable influence over their yield as 
well. However, the normalization of the spectra and thus particle yields 
are more directly controlled by the particle fugacity Tj = e"^'/^, where 
(Tj is particle 'i' chemical potential. Since for each related particle and 
antiparticle pair we need two chemical potentials, it has become convenient 
to choose parameters such that we can control the difference, and sum of 
these separately. For example for nucleons and antinucleons A^, N the two 
chemical factors are chosen as: 

CJAT = /ifc + rin7Ar, cr^ = -/ib + Tln7Ar, (3) 

TN=lNe^''^, T^ = 7^e-'^'/^. (4) 

The role of the two factors can be understood considering the first law 
of thermodynamics: 

dE + PdV-TdS = a^dN + ajjdN, (5) 
= UbidN - dN) +T In -fN{dN + dN). 

The (baryo)chemical potential fib, controls the baryon number, arising from 
the particle difference. 7Ar, the phase space occupancy, regulates the number 
of nucleon-antinucleon pairs present. 

There are many different hadrons, and in principle, we should assign 
to each a chemical potential and then look for chemical reactions which 
relate these chemical potentials. However, more direct way to accomplish 
the same objective consists in characterizing each particle by the valance 
quark content forming a product of chemical factors, e.g., for p{uud), 

Considering the isospin symmetry of strong interactions aside of the 
three quark fugacities = u,d,s, we also introduce the light quark fu- 
gacity: 

Ag = A„Arf, Xb = Xq. (6) 
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Table 1. Four quarks s, s, q, q require four chemical parameters; right: name of 
the associated chemical equilibrium 



7j controls overall abundance Absolute chemical 

of quark {i = q, s) pairs equilibrium 

Aj controls difference between Relative chemical 

quarks and antiquarks {i = q, s) equilibrium 



To transcribe the fugacities into chemical potentials we recall: 

Ai = e^»/^, ^q = ^{^u + l^d)- (7) 

The relation between quark based chemical potentials and the two principal 
hadron based chemical potentials of baryon number and hadron strangeness 
/ij, i = b,S is: 

/xft = 3/ig /^s = -^/^fe - /^s, Xs = ^- (8) 
6 As 

An important (historical) anomaly is the negative S-strangeness in s- 
hadrons, e.g.: 

There are two types of chemical factors ji and fii and thus two types of 
chemical equilibria. These are shown in tabled The absolute equilibrium 
is reached when the phase space occupancy approaches unity, 7^ — > 1. The 
distribution of flavor (strangeness) among many hadrons is governed by the 
relative chemical equilibrium. 

There is considerable difference in the dynamics of these two particle 
yield equilibration. This can be, e.g., understood considering strangeness 
in the hadronic gas phase. The two principal chemical processes are seen 
in Fig.^ The redistribution of strangeness among (in this example) A, tt 
and N, K seen on left in Fig.^ constitutes approach to the relative chemi- 
cal equilibrium of these species. The s, s pair production process, on right 
in Fig.^ is responsible for absolute chemical equilibrium of strangeness. 
Achievement of the absolute equilibrium, 7^1, require more rarely occur- 
ring truly inelastic collisions with creation of new particle pairs. 
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Fig. 1. Typical strangeness exchange (left) and production (right) reactions in the 
hadronic gas phase. 

2.3. Ratios of particle yields 

When particles are produced in hadronization, we speak of chemical 
freeze-out. To obtain the particle yields up to overall normalization the 
individual Fermi-Dirac and Bose-Einstein phase space integrals are evalu- 
ated. In order to arrive at the full yield, one has to be sure to include all 
the hadronic resonance decays feeding into the yield considered, e.g., the 
decay K* K + tt feeds into K and vr yields. This actually constitutes a 
book keeping challenge in study of particle multiplicities, since decays are 
contributing at the 50% level to practically all particle yields. Sometimes 
the decay contribution can be dominant, as is generally the case for the pion 
yield. For pions in particular, each resonance contributes relatively little in 
the final count, it is the large number of resonances that contribute which 
competes with the direct pion yield. 

To determine the parameters governing the chemical freeze-out, we an- 
alyze particle yields in terms of the chemical parameters and the temper- 
ature. Except for direct pions practically always one can use Boltzmann 
approximation and large reaction volume, and what follows in this subsec- 
tion assumes that this simple situation applies. 

It is often appropriate to study ratios of particle yields as these can be 
chosen such that certain physical features can be isolated. For example, just 
the two ratios, 

lead to a very good estimate of the baryochemical potential and strange 
chemical potential |^, and thus to predictions of other particle ratios. 

The sensitivity to phase space occupancy factors ji derives from com- 
parison of hadron yields with differing q, s quark content, e.g.: 

E-{dss) -fdls gsAdAg E'jdss) jd7s 9s^d^K^ 

A{dds) ^ -fjjs gA^lXs ' A{dds) ^ 7j7s 5AArf^A7^ ' 
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In Eq. , each of the ratios also contain chemical potential factors 
Aj. These can be ehminated by taking the product of particle ratio with 
antiparticle ratio, thus 



A{uds) A{dds) " V7«5A 

The proportionality constant describes the phase space size ratio for 
the two particles of different mass. It incorporates the contributions from 
resonance decays, which of course differ from particle to particle. 

The method applied in Eq. (|12() can be used in several other such double 
particle ratios. The relevance of this is that we have identified an experi- 
mental observable (combination of particle ratios) solely dependent on two 
parameters of statistical hadronization and chemical freeze-out, the temper- 
ature T which controls the phase space factor ratio C and the ratio "fs/lq- 
Another double ratio of considerable importance is 

A{uds) A{uds) _ ^2 fj±9Ay 

Both ratios Eqs. (|12I13|) are independent of chemical potentials, their mea- 
surement allows to constrain the value of 'js/^Jq only as function of T |18j . 

Double ratio made of mesons in general will be weakly dependent on 
chemical potentials since some kaons and pions are decay products of bary- 
onic resonances. The leading term is in general originating in direct and 
resonance mesons, for example: 

K^(u-s)K-(us) clJ:^i^)\.... (14) 



TT+{ud) p (ud) \-fg 

In summary, and for the benefit of our later discussion in subsection 
15.11 we remember that there is no dependence in Eqs. H12I13|1 and almost 
none in Eq. H14|) on the baryochemical potential and strange chemical po- 
tential. This eliminates dependence of the baryon and strangeness density, 
and allows to focus on the physics issues unrelated to baryon compression 
effects. 

2.4- Particle phase space density and hadronization 
The maximization of microcanonical entropy, 

Sf,B = I^^T[{lTfF,B)HlTfF,B)-fF,BlnfF,B] , (15) 



10 



zakopOSc printed on February 1, 2008 



(where indeed minus sign is for fermions and plus sign for bosons) subject 
to energy and baryon number conservation implies the quantum distribu- 
tions pij : 

(f'Ni Qi 1 ^bosons ^ ^rUi/T ^-j^g-j 



When the phase space is not densely occupied, ±1 in the denominator can 
be neglected and we have the Boltzmann approximation we used in last 
subsection: 

^ -^i —J7;I'J^ ('\'7\ 

dW^"^'(2^^ ' • ^ ^ 

In Eq. (fT6|) . for the Boson distribution, the value of 7g is limited by the 
condensation singularity. The maximum value of ^^'^^ = g»"'r/(2T) pj^ys a 
very pivotal role considering that the mass of the pion and hadronization 
temperature are similar. Large value of 7g — > e™""/"^ can be directly noticed 
in pion spectra in an uptilt in the soft portion of the distribution. A 
similar constraint is also arising from kaon condensation for 7^ but it is 
much less restrictive: 



7g V-^- 



±1 

< e^^K-™-^/' ~ 11 . (18) 



In the local restframe the particle yields are proportional to the momen- 
tum integrals of the distribution Eq. (jlGj) . As example, for pions vr, nucleons 
A*" and antinucleons we have: 

N = CVgNf 7^ (21) 

^%(27r)3i + ^-3A+3ei?/T ^ ' 

Naively, we would think that the coefficient is simply the volume V, as it 
would be the case for a gas of hadrons. However, more generally there is 
a common additional factor C which is determined by the dynamics of the 
hadronization process. The region of integration R is determined in terms of 
the experimental acceptance, keeping in mind that we are here considering 
the phase space in rest frame of hadronizing QGP, while the experimental 



zakopOSc printed on February 1, 2008 



11 



detector is at rest in the laboratory. Both the motion of the cohision center 
of momentum frame with respect to the detector, and the collective flow of 
QGP have to be allowed for. 

Generally a small region of particle rapidity y but nearly all of range of 
is accepted in experiments at RHIC. Assuming the yield of particles is 
practically constant as function of y, it is possible to imagine that the yield 
arise from a series of fireballs placed at different rapidities. Thus in this 
limit we can proceed to study hadronization as if we had a full phase space 
coverage. However, in this case in particular the proportionality constant 
C can be quite different from unity. 

The hadronization of QGP phase is not fully understood. However, 
numerous experimental indicators show that hadrons are emerging rapidly 
from a relatively small space-time volume domain. In such hadron 
formation has to absorb the high entropy content of QGP which originates in 
broken color bonds. The lightest hadron is pion and most entropy per energy 
is consumed in hadronization by producing these particles abundantly. We 
evaluate Eq. ()19() to find particle number, and the entropy content follows 
from Eq. (fTC)): 

[(1 + /.) ln(l + - U In U] , (22) 

UE) = _2 J-,^ ^ , E^ = ./^^;^T^. (23) 

7g e^-^/-^ — 1 

As is seen in Fig.|^ the maximum entropy density S/V occurs for an 
oversaturated pion gas, 7q ~ Qm^r/^T ^ -j^ g_ Here, the entropy density of such 
a saturated Bose gas is twice as large as that of chemically equilibrated Bose 
gas. Since aside of pions also many other hadrons are produced, the large 
value of 7g is necessary and sufficient to allow for the smooth in iJ,b,T,V 
transformation of a QGP into hadrons. The number of active degrees of 
freedom in the oversaturated hadron gas with jg — > and in 'freezing' 
QGP phase is very similar. 

3. Strangeness — a popular QGP diagnostic tool 

3.1. Motivation for study of strangeness 

It is widely agreed that the process of strange particle production in rel- 
ativistic nuclear collisions offers considerable insight into the structure and 
dynamics of the dense matter formed in these reactions. The experimental 
reasons for the interest are practical, as this observable is experimentally 
accessible. The theoretical interest has been motivated by the recognition 
that strangeness is a marker of the thermal gluonic degree of freedom of the 
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deconfined phase [13 HH] • These result is enhanced abundance of strangen- 
ess, and thus production of numerous strange hadrons j29j . 

In an nutshell the experimental motivation for study of strangeness is: 

• There are many strange particles allowing to study different physics 
questions {q = n, d): 

(piss), K{qs), K{qs), A{qqs), A{qqs), R{qss), E{qss), Q,{sss), Q,{sss); 

• Strange hadrons are subject to a self analyzing decay within a few cm 
from the point of production. As an example, we see the cascading 
decay of a doubly strange H~ in Fig.El 

• Though some of the strange hadrons are produced quite rarely in p-p 
reactions, the enhancement we encounter in A- A reactions produces 
'exotic' strange observables with yields sufficient to reach a relatively 
high statistical significance. 

The theoretical case for study of strangeness has been researched in 
depth: 



zakopOSc printed on February 1, 2008 



13 




Fig. 3. S -decay, dashed line the invisible A emerging from the decay 
kink ending in the decay 'V of the final state charged particles. 



Fig. 4. Lowest order Feynman diagrams for the production of flavor. Top line: 
gluon fusion and bottom line: quark-antiquark fusion into strange quark pair. 

• Production of strangeness occurs predominantly in gluon fusion gg — > 
ss. Thus abundant strangeness is linked to presence of thermal gluons 
from QGP. The quark process 30, , — > ss, contributes at level of 10- 
15% of total rate, the Feynman diagrams for the lowest order process 
are seen in Fig.|lJ 

• Coincidence of scales: ~ implies that the relaxation time re- 
quired to equilibrate strangeness yield chemically, r^, is of the same 
magnitude as the lifespan of QGP, Tg ^ tqgp- In consequence, strang- 
eness acts clock for the precise QGP lifespan. 

• Due to presence of baryon density, the yield of light antiquarks is sup- 
pressed compared to yield of quarks and one thus can encounter s > q. 
Ensuing recombination production of antibaryons at hadronization of 
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QGP can lead to a relative strange antibaryon enhancement, and at 
RHIC to (anti)hyperon dominance of (anti)baryons j31j . 

3.2. Two step hadron formation mechanism in QGP 

The establishment of a ready supply of strange quarks occurs in a man- 
ner independent of the production of final state hadrons. After some time 
has passed, QGP drop reaches hadronization condition in its internal pres- 
sure driven expansion. Then, the final state hadrons emerge, as is seen in 
Fig-El Thus strange hadron production occurs in two independent steps 
following each other in time: 

1. gg ^ ss predominantly in the early hot QGP; 

2. hadronization of pre-formed s, s quarks after QGP cools. 



Fig. 5. Illustration of the cross-talk two step mechanism of strange hadron for- 
mation from QGP: inserts show gluon fusion into strangeness, followed by QGP 
recombinant hadronization. 

This sequence allows the production of complex rarely produced multi 
strange (anti) particles, enabled by 'cross talk' between s, s quarks made in 
different microscopic reactions. We consider observation of this process to 
be a specific signature of deconfinement, since it provides evidence for both 
the activity of thermal gluons, and the mobility of strange quarks. The 
specific feature of such an enhancement is that the enhanced production of 
strange antibaryons is increasing with strangeness content. 
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3.3. Strange hadron enhancement or suppression? 

The enhancement of strange hadrons has been studied as function of the 
number of participating baryons, which are evaluated in terms of 'wounded 
nucleons', that is nucleons which have been involved in an inelastic scat- 
tering in the collision. We show, in Fig. El the results of NA57 |32| I33j . 
which extend those presented by WA97 collaboration. The yield of strange 
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Fig. 6. Specific yield per wounded nucleon of strange baryons A, S and antibaryons 
A, S and of 17 + as function of the number of wounded nucleons [SSI- The data 
are obtained in 5 centrality bins which cover the range 50-350 participants. The 
baseline for the enhancement is the yield observed in the p-Be reactions. 

hadrons shows the expected increase of enhancement as function of newly 
made quark content s, s, q, and as function of the participant number (see 
for example the result shown in figure 37, in j34j . The enhancement is in 
Fig. reported with a base for p-Be reactions. It has been recognized long 
ago that the relatively small number of strange quarks produced in proton 
induced reaction could require evaluation of statistical yields applying mi- 
croscopic particle number constraints of relevance in such a limit • Thus 
statistical equilibrium yields for the p-p reaction systems are suppressed and 
thus enhancement is seen comparing AA to pp. However, results of NA57 
we see in Fig.|Hl disagree with the expected behavior of this suppression pS] . 

Without going into the theoretical details discussed elsewhere the 
important feature of the micro-canonical yield suppression is that the asymp- 
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totic particle yields are relatively quickly attained as is seen in Fig.d We 
see here the enhancement Ei generated by rebasing the small system yields 
to unity for the strangeness content i = 1,2, 3. The top of the figure shows 
the resulting enhancement for the p-p reaction system while the bottom 
estimates the result for p-Be as used by NA57, assuming that the yield of 
strange quarks pairs {Ng) doubles from 0.66 to 1.3. The size dependence is 
expressed in terms of the canonical yield of strange quarks produced (thus 
the base enhancement unity is at the value A'^^^ = 0.66 in the top portion 
of Fig.[7| while it is at N'^^ = 1.3 in the bottom portion. 

102 
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Fig. 7. Enhancement Ei of strangeness i — 1 (solid lines), i = 2 (long dashed 
lines) and i — 3 (short dashed lines) hadrons as function of the canonical yield of 
strangeness N'~^^ . The cases for p-p (top panel) and p-Be (bottom panel) reactions 
are shown. 

Looking at the results in Fig. El we see that systems with 100-300 
baryons are still growing in yield. On the other hand we see in Fig. [3 that 
the saturated yield of strange hadrons is attained at A^*^^ = 5, which yield 
is attained for a fireball with 10-30 baryons. This, along with the other 
issues |37| . raised earlier shows that the enhancement effect is indeed not 
due to a suppression of the base of comparison, but is due to kinetic pro- 
cesses already described, i.e., excess strangeness production in QGP and 
subsequent hadronization process. 
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3.4- Sudden hadronization 

A fireball is not a piece of deconfined matter sitting still. Given its origin 
in a collision on two atomic nuclei, it is initially made of highly compressed 
matter. Thus it subsequently undergoes a rapid collective explosive outflow. 
Such a collective motion of color charged quarks and gluons contributes 
an important collective component in the pressure, beyond the rest frame 
thermal pressure. This can be seen considering the stress portion of the 
energy-momentum tensor: 

=P5,, + {P + e)r^^. (24) 

The rate of momentum flow vector V at the surface of the fireball is obtained 
from the energy-stress tensor T^i: 

p = f-n = Pn + {P + e)^^. (25) 

The pressure and energy comprise both the particle and the vacuum prop- 
erties: 

P = Pp-B e = ep + B. 
The condition V = reads: 

Sn = Ppn + (Pp + ep)^^. (26) 

Multiplying with n , we find: 

^ = Pp + (Pp + ep)-^, K=^^^. (27) 

We note that the hadronization hypersurface, and in particular the angular 
relation between the collective flow direction vector and surface normal 
direction, defines the force balance condition seen in Eq. ()27() . 

In order to satisfy the condition Eq. (|27|) . we must have Pp < B: the 
QGP phase pressure P = Pp — B must be negative at hadronization since 
it has to compensate the effect of the positive dynamical pressure. When 
P — > and a nonzero velocity of collective flow remains pointing outward, 
the surface is torn apart in a rapid fllamentation instability. This situation 
can only arise since the quark-gluon matter presses again the collective 
vacuum which is not subject to collective dynamics [38j . In this aspect the 
dynamics of QGP expansion reminds us of a gas bubble growth within a 
liquid (cavitation phenomenon). 
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The hypothesis that hadronization of the QGP deconfined phase formed 
in high energy nuclear colhsion is sudden has a relatively long history. In 
the 1986 review Koch et al. [SHI note that strange hadrons will be viable 
observables of QGP should a long lived hadron reequilibration not occur. 
They offered a study of hadron abundances arising from sudden hadroniza- 
tion of strangeness rich QGP. Once data became available, the first analysis 
of the experimental hadron abundances is carried out within this approach 
j26j . The important feature of sudden hadronization, the identity of tem- 
peratures of diverse hadrons and in particular of baryons and antibaryons is 
recognized. Shortly after dynamical models of sudden hadronization begin 
to emerge |^ • 

Phase boundary for a system at rest between the hadron gas domain 
and quark-gluon plasma in the Hb,T plane, is shown in Fig.|Hl The solid 
thin line is that for hypothetical point hadrons, while dashed line shows the 
typical estimate for finite volume hadrons. The effect of the 'wind' of fiow 
of QCD matter is also seen in Fig. El for a geometry of hadronization with 
K = 0.6. 



1000 

800 

^ 600 
3. 

400 
200 


100 110 120 130 140 150 160 170 180 
T [MeV] 

Fig. 8. Hadronization boundary in the fi^, T plane, QGP-HG transition for point 
hadrons (solid thin line) and finite volume hadrons (dashed line), are typical es- 
timates of QGP-HG boundary for a system at rest. Dotted line for finite size 
hadrons includes the effect of the collective flow velocity for Vc — 0.32, k — 0.6. 
Thick solid line: breakup with v — 0.54, k — 0.6 (based on |38| ) . 

As the expansion velocity increases, the magnitude of the color wind 
effect increases with velocity and is here shown for Vc = 0.32 (dotted) and 
Vc = 0.54 (solid lines). With increasing fiow velocity the phase boundary 
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moves to lower temperatures. This effect reduces for a given ^Uf, the magni- 
tude of critical temperature in a significant way. 

As collective flow velocity decreases, the hadronization temperature in- 
creases. Even though this effect is somewhat compensated by the influence 
of the increase in ni, which implies a decrease in T, this 'normal' curvature 
effect seen in each of the lines in Fig.|Hlis smaller than the 'new' flow effect as 
we move from line to line in Fig.|Hl At RHIC, the largest transverse velocity 
is reached, and thus largest supercooling achieved. We expect the largest ef- 
fect from the wind of QGP flow, leading to smallest observed hadronization 
temperature. 



We use the statistical hadronization model presented in section [21 to fit 
experimental particle yields and to evaluate the hadronization conditions 
including the specific strangeness, energy and entropy content. Experimen- 
tal results are now available for several collision energies both at SPS and 
RHIC. Thus it is possible to study the energy dependence systematically 
j42j . We will see if there is consistency in the fitted energy dependence 
both within the energy range of SPS and RHIC, and between these two dif- 
ferent energy domains. Specifically, we reevaluate here our comprehensive 
analysis of RHIC-130 Au-Au colhsions [HIISI, carry out a RHIC-200 pre- 
liminary analysis and compare with SPS Pb-Pb stationary target reactions 



at \ s%^^ = 8.75, 12.25, 17.2 GeV (projectile energy 40, 80, 158A GeV). 



For SPS, we use here the CERN-SPS-NA49 47r particle multiplicity re- 
sults HmS], which include 7r±, K, K, A, A, (p at 40, 80, 158^ GeV. We 
also fit (relative) yields of H, H, 0, when available. Our current results 
are found quite near to our earlier work reported for the 158^ GeV energy 
j43j . even though our data sample has been limited to (a subset of) NA49 
47r-results, while our earlier work was dominated by WA97 central rapid- 
ity results. Since we fit only 47r-particle yield, no information about the 
collective flow velocity is obtained. 

Our statistical hadronization program allows for decays of resonances, 
and for non-equilibrium occupancies of light quarks "jq ^ 1 and strange 
quarks 7<j 7^ 1. In the results, we discuss here, we enforce strangeness con- 
servation. We assume that 50% of weak decays from H to A and from to H 
are inadvertently included in the yields when these had not been corrected 
for such decays. Particle yields containing this type of contamination from 
a weak interaction cascade are marked by a subscript 'c'. We assume that 
the pions from such weak decays are not included in the experimental yields 



4. Hadronization at SPS and at RHIC 



4.1. Method of analysis 
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as these pions tracks clearly do not originate in the interaction vertex. 

Not all pion producing resonance decays are known, and many have to 
be guessed given the present status of the particle data tables. We have 
estimated comparing to expectations based on an exponentially rising mass 
spectrum that pion yields may be undercounted in the data based statistical 
hadronization programs |13j . However, one can also argue that pion yield 
is overcounted: given the sudden nature of the hadron fireball breakup, 
microscopic reactions cannot exactly follow through to populate the heavy 
resonances. To estimate the potential for systematic error we study the 
hadronization with pion yield artificially changed by ±15%. We find, not 
surprisingly, that there is an anti-correlation to fitted temperature in the 
range of =f5%, respectively, but other statistical parameters and physical 
properties remain stable. 

Some of the particle ratios we consider can be formed from other results 
already considered. However, each measurement has to be considered as 
being independent and often the data which theoretically is related by a 
product of ratios, is in praxis barely consistent. The usual recommended 
procedure which we follow here is to include all the available different par- 
ticle ratios in the fit, and to include for each measured data point its error 
in the study of statistical significance. When several experimental results 
are available for the same ratio we fit the average result. 

4-2. RHIC analysis 

In the study of particle yields at RHIC, we have allowed for both sta- 
tistical and systematic errors in the experimental data set. We show the 
results of the fit procedure at RHIC in table|21 The top line gives the re- 
action energy. The fitted statistical parameters follow, with errors being 
both statistical and systematic, derived from variation of the pion yield. 
The equilibrium results for statistical parameters and particle multiplicities 
agree very well with those presented for the case of chemical equilibrium by 
the Krakow group [45] . 

As the bottom section in table[21shows, the statistical significance is very 
good for chemical non-equilibrium fit, and also, we note that total is half 
as large as it is for the equilibrium case on the right. Only for a very large 
number of degrees of freedom one converges toward x^/dof=l when the 
experimental errors are correctly evaluated. For a relatively small number 
of degrees of freedom, as is the case in particular for RHIC-200, we must have 
X^/dof significantly smaller than unity to believe in the physical significance 
of the result. The bottom line presents an estimate of the probability Ptrue 
using the values of and dof that the fit is an appropriate physics model 
description. Since the systematic errors are here taken into account, this is 
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Table 2. The RHIC chemical freeze-out statistical parameters found for non- 
equilibrium (left) and semi-equilibrium (right) fits to RHIC results. We show 
^Jsnn, the temperature T, baryochemical potential /i;,, strange quark chemical 
potential /i^, strangeness chemical potential /is, the quark occupancy parameters 
7q and js/lq, and in the bottom line the statistical significance of the fit. The star 
(*) indicates that there is an upper limit on the value of 7^ < e™"/"^ (on left), 
and/or that the value is set (on right). 



^SNN [<^eV] 


200 


130 


200 


130 


T [MeV] 


143 ±7 


144 ±3 


160 + 8 


160 ±4 


/^b [MeV] 


21.5 ±31 


29.2 ±4.5 


24.5 ± 3 


31.4 ±4.5 


f^s [MeV] 


2.5 ±0.2 


3.1 ±0.2 


2.9 ±0.2 


3.6 ±0.2 


/xs [MeV] 


4.7 ±0.4 


6.6 ±0.4 


5.3 ±0.4 


6.9 ±0.4 




1.6 ±0.3* 


1.6 ±0.2* 


1* 


1* 




1.2 ±0.15 


1.3 ±0.1 


1.0 ±0.1 


1.13 ±0.06 


xVdof 


2.9/6 


15.8/24 


4.5/7 


32.2/25 


-ftrue 


90%+ 


95%+ 


65% 


15% 



an appropriate time to look at this value. Clearly as more data for RHIC- 
200 becomes available, and errors shrink, this evaluation wih change. 

One is often interested to see how well the particle yields fit the exper- 
iment. For RHIC-130, we obtain the fitted yield ratios of the following 24 
results, along with a few (three) results one would wish to also have avail- 
able (see table 131) • It is important that the reader realizes that these ratios 
can depend significantly on the assumptions made about the cascade of 
particles that ensues hadronization, and which we cannot fully characterize 
given lack of experimental data. We have 'played' around with branching 
ratios of high mass resonances and believe that the changes that ensue are 
all within the errors given for the statistical parameters. However, the indi- 
vidual yields presented below may change, and thus other groups can derive 
slightly different ratios for same statistical set of parameters. 

For RHIC-200 we employ the following 10 relative yields available at this 
time (see table We see that the presented non-equilibrium multiplicity 
fits are good, and the individual contributions to are satisfactory. 

The quality of RHIC results presented in table[21confirms that the statis- 
tical hadronization mechanism is an appropriate particle production model 
in presence of the sudden QGP breakup which without doubt is seen at 
RHIC. However, the statistical hadronization approach may be not entirely 
appropriate at low SPS energies, and all AGS energies. It is possible that 
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Table 3. 



ratio 


RHIC-130 


non-eq fit 


x' 

A, 


TT IP 


ij .O _1_ 1. . '-t 


Q 07 


1 1 5 


" If 


13.4 ± 0.9 


13.15 


0.08 






0.0459 






059 ± 004 


0509 


4.11 


J. j 1 V 


0.042 ± 0.004 


0.0379 


1.04 


c / 


0.0079 ± 0.0012 


0.00805 


0.01 


-~ Ih- 


n nnfifi + n nni 

u.uuuu in u.uux 


n nnfi4'i 

U.UUU'itJ 


n 09 


o //)- 


n 9 -1- '^'ii 


1 9 1 


u.uu 


\\l -\- \ L) / n 


\ll ± D.O jiU 


z4.o iU 


U.iy 




n Qn 4-019 
u.yu ± u. iz 


7/17 


i.Do 




U.yo it u.iy 


U.OZD 


U.oU 




n 1 no _i_ n no 
O.i9o ± U.Uo 


0.189 


n no 
O.Oz 


1 — ' /A 


0.219 ± 0.035 


0.207 


0.12 


i £ / 1 1 




n 1 R4 








0.180 




(n+f])/(H + H) 


0.150 ±0.04 


0.171 


0.28 


p/p 


0.71 ±0.06 


0.674 


0.36 


Ac/Ac 


0.71 ±0.04 


0.745 


0.78 




0.83 ±0.08 


0.801 


0.13 




0.95 ±0.1 


0.878 


0.51 




0.17 ±0.02 


0.195 


1.59 




0.15 ±0.02 


0.180 


2.28 




0.87 ±0.07 


0.923 


0.57 




0.26 ±0.08 


0.231 


0.13 


ct^/h- 


0.02 ± 0.002 


0.0212 


0.37 




0.15 ±0.03 


0.148 


0.00 


(t)/K*^ 


0.595 ± 0.24 


0.639 


0.03 



kinetic models need to be applied, which ahow for a multitude of freeze-out 
conditions depending on the nature of particle considered. For this reason, 
when we study this energy domain we will look more closely at the behav- 
ior of individual observable {K/ir) rather than at the global yields of all 
particles, the method we are using at RHIC and high SPS energies. 

4-3. Phase space occupancy, energy stopping as function of collision energy 

The non-equilibrium parameters Js^lq we find in our fits are shown in 
Fig.|^ as function of collision energy. In top figure, open squares show the 
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Table 4. 



ratio 


RHIC-200 


non-eq fit 


x' 

A. 


TT IP 




Q 8fi 




TT jp 


1 9 1^ -1- 1 7 


lO.U 




vlhr 




04759 




■t-^Cj til 




0.0441 




7C/h~ 




0.0354 




rr Ihr 




00657 








U.UUooo 




t u — 

il/h 




iU ID 


0.31 






U.Do / 




A / - 




0.738 




" /A 




0.175 




/A 




U.io ( 








0.152 








0.162 




p/p 


0.74 ± 0.04 


0.746 


0.02 


Ac/Ac 




0.801 








0.844 






1.05 ±0.2 


0.902 


0.55 






0.182 




K- jix- 


0.156 ±0.02 


0.172 


0.63 


K-/K+ 


0.95 ± 0.05 


0.945 


0.01 


K*^/K- 


0.205 ± 0.033 


0.227 


0.44 


4>/h- 


0.02 ± 0.002 


0.0183 


0.74 


<t>lK- 


0.118 ±0.04 


0.133 


0.14 




0.595 ± 0.123 


0.585 


0.01 



result when we fix 7^ = 1. We note that in this case the SPS fits yield 
7s 2± 0.7. However, we see that 7^ > 1 is preferred in all cases, with value 
converging toward the limit 7^ — > 7™^ = g^i^/sT 'pj^^jg behavior reflects 
on the need to account for the excess in charged hadron multiplicity. With 
this large 7^ also 7s is found to be quite large, in particular so at RHIC. We 
need to solidify this result since this is indeed the expected signature of the 
sudden hadronization of deconfined strangeness rich QGP phase. 

We next address the reaction mechanisms, evaluating the energy stop- 
ping. We define it as the ratio of the thermal energy Ef^j^ found in particles 
produced and normalized to a nucleon pair, with the collision energy -^snn 
which is also given per colliding nucleon pair. This study will allow us to 
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Fig. 9. Fits of, top to bottom 7s (top) and 7g (bottom) chemical non-equilibrium 
parameters as function of collision energy. Open squares results for 7^ = 1. 



consider the energy dependence of other variables as function of the energy 
available E^^j^ in the fireball, rather than as function of the energy of col- 
liding nuclei y^SAW; which in its greatest part may not participate in the 
reaction. 

The open triangles Fig. represent energy stopping for the fit assuming 
full chemical equilibrium. We note that at SPS this leads to the absurd be- 
havior that stopping is increasing with increasing collision energy, reaching 
70% at the top energy. This alone proves that the chemical equilibrium 
approach to statistical hadronization is physically meaningless. 

The open squares in Fig.^J show the stopping found for the case of 
semi-equilibrium fit, i.e., allowing only strange quarks to be out of chemical 
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Fig. 10. Fraction of energy stopping at SPS and RHIC: results are shown for 40, 
80 and ISSA GeV Pb-Pb fixed target SPS reactions and for 130 and 200A GeV 
Au-Au RHIC interactions. 



equilibrium. The full squares in Fig. apply full chemical non-equilibrium 
with phase space occupancy shown in Fig.l^l Up to the step-up between SPS 
and RHIC we see a smooth fall of of the stopping with increasing collision 
energy when chemical non-equilibrium is allowed for. We recall that the 
statistical hadronization study we have made assumes tacitly that Bjorken 
rapidity scaling applies exactly. One would hope that the RHIC stopping 
result will come down when we allow for the non-scaling yields of particles. 



4- 4- Strangeness and entropy in statistical hadronization 

We would like to compare the efficiency of strangeness and entropy pro- 
duction between the two widely different systems considered. We evaluate 
the number of strange quark pairs produced (strangeness yield) and divide it 
by the computed thermal fireball baryon number. This eliminates the need 
for an absolute yield normalization. More generally, a ratio of two extensive 
computed variables is nearly independent of the dynamics of hadronization. 

In Fig. 1111 we present strangeness production per baryon. The top 
Fig. ^2 section shows our result as a function of the collision energy. When 
using the thermal energy content in the reaction, we arrive at the result 
with a smoothly rising curve seen in the bottom figure section. As noted in 
the insert, strangeness yield is rising faster than linearly with fireball energy 
content. 
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Fig. 11. Strangeness per thermal baryon participant, s/h as function of ^/sJjn 
(top), and as function of thermal specific energy content E^^j^^ (bottom). 



We next look at strangeness per entropy, s/S, shown in Fig. 1121 There is 
considerable physical importance of s/S, since both entropy S and strang- 
eness s are produced early on in kinetic processes, are nearly conserved in 
the hydrodynamical expansion of QGP, and increase only moderately in 
the hadronization of QGP. Thus the value seen for s/S is reflecting on the 
kinetic mechanisms operational in early stages of the heavy ion collision. 
Because there is strangeness mass as energy threshold to overcome, one 
would expect that as the initial QGP phase becomes hotter with increas- 
ing collision energy, the excitation of strangeness grows more rapidly than 
excitation of entropy. 

We note, in Fig. 1121 for the chemical non-equilibrium assumption, a 
smooth and slow increase of s/S, from 0.02 to 0.025 at SPS. The RHIC- 
130 point is within the smooth extrapolation of the SPS result. The low 
value of RHIC-200 result is probably expressing a systematic error related 
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Fig. 12. Strangeness per entropy s/S as function of iJfJvAf- 



to the absence of multistrange particles in the fit. While entropy is already 
fully accounted for, only availability of multistrange hadron results will al- 
low to obtain a reliable strangeness yield. However, if this expectation is 
not fulfilled, it will be very important to obtain experimental results in the 
energy range between SPS and RHIC, in the domain where the change of 
the systematic behavior (slope) of s/5 as function of energy occurs. 



5. Search for an Energy Threshold 

5.1. Kaon to pion ratio 

One of the most interesting questions is if there is an energy threshold 
for the formation of a new state of matter. In our study of strangeness per 
entropy s/S, in subsection 14.41 we have seen a possibility that new physics 
sets in between the SPS and RHIC-130 energy domain. However, a final 
verdict on this issue depends on further RHIC-200 results. 

The NA49 experiment has been exploring the low energy domain at SPS 
in search for a new physics threshold between AGS and SPS energy range. 
We refer to the report by Marek Gazdzicki for the set of arguments that just 
at the bottom of SPS energy range something new happens [ITj. Perhaps 
the most important evidence cited by Marek is a peak in the /tt^ ratio, 
see top section of Fig.^J [15| . 

We think that this peak is purely an effect of the dense baryon medium. 
As the collision energy increases, beginning with low energies, the baryon 
density increases. However, at some collision energy baryons shoot through 
and the baryon density begins to drop. Since there is considerable sensitivity 



28 



zakopOSc printed on February 1, 2008 




0,2 



0.1 



+ 



B 



v's-lROO- 



10 100 
[GeV] 

Fig. 13. Top: 47r-ratio of yields of /n^ for nuclear (filled symbols) and ele- 
mentary interactions (open symbols); (courtesy of NA49 collaboration 45, 52]). 
Bottom: the reduced K/tt ratio for nuclear (filled) and elementary collisions (open 
symbols) both as function of collision energy. Figures rescaled to the same magni- 
tude on abscissa and ordinate. 



to baryon chemical potential and thus baryon density in the /tt^ ratio 
considered by the NA49 collaboration, we prefer to consider the nearly 
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baryon density independent ratio Eq. (|14jl , 



K K+ K- , , 

- = \/— — • (28) 

vr V 'K^ TT 

We present K/tt double ratio in the bottom of Fig. as black filled 
square. At high energy, we show the two RHIC results and at SPS and 
AGS energies we use the NA49 data set 48, 52j. The open squares are 
for the charged ratio K+/7r+ measured in pp reactions, which is offering 
an upper limit on K/tt from pp reactions, also using NA49 data set. We 
indicate for y/s = 1800 the p~p TEVATRON result ^51.. 

Both upper and lower portion of Fig. are drawn on same scale. Com- 
paring the top and bottom in Fig.^1 we see that the peak between AGS 
and SPS energy domain has completely disappeared, and that in the SPS 
energy domain the K/tt ratio Eq. (|28jl is flat. This behavior we already saw 
in the ratio of strangeness to entropy, bottom of Fig. 1111 Moreover, when 
we look at the inverse slopes of the K"*" mj_-distributions. Fig. 1141 we again 
see in SPS domain a flat distribution [15] . 
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Fig. 14. The inverse slope T of K+ spectra at AGS (filled triangles) at SPS 
(filled squares, including a recent 30A GeV NA49 experiment and at RHlC-130 
and RHIC-200. 

Although the large error bars leave a lot of space for further structure, 
we do not see any evidence for a 'peak' in strangeness production, rather 
we believe that NA49 has discovered that K^ /tt^ ratio allows to determine 
where the baryon density peaks as function of collision energy. 
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As the dashed hnes in the bottom portion of Fig. El suggest, there is a 
smooth transition from rising to saturated behavior of K/tt ratio. However 
this behavior mirrors the result we see in the pp reactions, except that it 
occurs at lower collision energy. We conclude that if kaon and pion yields 
are solely used as the measuring stick for new physics, there is clear evi- 
dence that new physics is seen at RHIC, where a true deviation from pp, pp 
behavior occurs. 

5.2. A 'conventional' K/ti excess 

The enhancement of K/tt ratio at RHIC may be, however, simply result 
of novel mechanism of pion suppression. We first recall that there is con- 
siderable increase in the transverse velocity of matter flow, as is born out, 
e.g., in the higher inverse slopes of particle spectra at RHIC. We see this 
effect in Fig.Ulfor K+ . 

The kinetic energy of this transverse motion must be taken from the 
thermal energy of the expanding matter, and ultimately this leads to lo- 
cal cooling and thus a reduction in the number of quarks and gluons. The 
local entropy density decreases, but the expansion assures that the total 
entropy still increases. Primarily, gluons are feeding the expansion dynam- 
ics, while strangeness being weaker coupled remain least influenced by this 
dynamics. Model calculations show that there is practically no strangeness 
reannihilation in transverse expansion at RHIC •").{ . 

The depletion of the non-strange degrees of freedom in the feeding of 
the expansion assures an increase in the K/t: ratio with increase of collision 
energy. This effect explains why 7^ > 7^, see table If depletion by 
transverse expansion dynamics of the non-strange hadrons is the origin of 
the increase in K/t: at RHIC, then this effect would arise gradually as 
function of collision energy. If, however, new phase of matter is the cause 
for this increase, we would expect a more 'edgy' onset of the increase in the 
energy domain between SPS and RHIC. 

An even greater effect can be expected at future CERN LHC collider, for 
heavy ions LHC will be a GOOOGeV collider, 30 times the maximum energy 
at RHIC. The transverse expansion dynamics effects should be much greater 
at LHC and we expect significantly grater K/tt ratio, and thus greater value 
of ^sMq- The only limit that we can see for a rise in 7s /7g, which is the 
factor controlling the rise in K/t:, is the condensation limit, Eq. H18() : 

smce < e T and js7q < e t we have Is/lq < e t 

We have shown the resulting maximal value of K/tt in Fig.^Jas a solid line 
as function of hadronization temperature. 
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Fig. 15. Kaon to pion ratio K/t: (see text) as function of hadronization temperature 
T for chemical equilibrium (dashed line) and with maximum allowable 7^ and 7^ 
(solid line). 

The dashed hne in Fig.H^l shows the result when 7s = 7g = 1. Depend- 
ing on what we beheve to be a vahd hadronization temperature for a fast 
transversely expanding fireball, the possible enhancement in the K/tt ratio 
will be in the range of a factor 2-3. The arrow indicates whereto the ratio 
converges for T ^ 00. This limit applies to both solid and dashed lines and 
critically depends on the high mass resonances which we really do not know 
very well. 

It can be assumed that hadronization at LHC occurs in a temperature 
domain T = 160 it 20 MeV. In this case even more interesting than the K/tt 
ratio enhancement would be the associated enhancement anomaly in strange 
(antibaryon) yields. One would find a strong inversion in the population, 
with the more strange baryons and antibaryons being more abundant than 
less strange species. 

5.3. Final remarks 

We have shown that strangeness at CERN-SPS and BNL-RHIC is a well 
developed tool allowing to study QGP-hadronization. Our discussion has 
shown that a systematic study of strange hadrons fingerprints the properties 
of a new state of matter. We have argued that the deconfinement specific 
observable is strangeness, and we have explained why we have little doubt 
that deconfined phase has been formed in relativistic heavy ion collisions. 
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The situation about deconfinement onset as function of energy is today 
much less clear than ever before. We argue here (subsection I5.1|) that be- 
tween AGS and SPS energy range a smooth transition to a saturated K/tt 
yield is seen. This behavior mirrors the behavior of slopes of kaon spec- 
tra. At SPS, we note a great enhancement of strange antibaryon yields, 
rising with strangeness content. Analysis of strange hadrons shows that 
the specific strangeness yield is very high. The strangeness production and 
hadronization at SPS is consistent with the behavior of the QGP phase. 

Somewhere between SPS and RHIC the yield of K/tt begins to rise above 
the elementary reaction result. We saw, in the same energy range, a possible 
anomaly in the strangeness to entropy ratio s/S, subsection 14.41 However, 
we argue in subsection 15.21 that the K/ vr ratio rise could be a dynamical 
effect due to rapid transverse expansion, and we discuss in subsection 14.41 
that the s/S anomaly could be removed by more comprehensive results 
forthcoming from RHIC-200 run. Thus, at this time, RHIC could harbor 
some great new discovery, or be a QGP producer with a more powerful 
transverse expansion and grater energy density. The behavior of particle 
spectra support this view 

The BNL-RHIC experimentalists announced recently the discovery of 
new physical phenomena, jet quenching, as would be expected in QGP. In 
February 2000 CERN-SPS experimental groups have combined their exper- 
imental results to claim a discovery of new state of matter, which looks 
like QGP. Aside of strange hadron evidence, several other observables were 
implicated in the CERN announcement. Discovery of QGP is also claimed 
in the study of pp reactions |55j . This work evaluates the ideal gas degrees 
of freedom which number is enhanced as expected in QGP phase. 

Though there was not much AGS public relation campaign, we do not 
see a reason why deconfined state formation did not occur at these low 
energies. If the criterion of relevance is the heating of nuclear matter to 
the phase boundary found in lattice gauge calculations, this is achieved 
at the mid-range of AGS collision energy. However, in this case the new 
phase of matter is embedded in the center of a opaque hadron shell. Only 
electromagnetic probes would have allowed to see into the heart of such low 
energy fireball. A new experimental facility with much more intense beams 
is needed to perform such experiments. 

A student reading this report may find a useful summary of facts about 
QGP, physics of strangeness, and statistical hadronization. The advanced 
QGP researcher will find here practical and useful results addressing analysis 
of the RHIC-130 and RHIC-200 experimental data, and a comparison with 
SPS results. These results do suggest that through study of strangeness at 
RHIC we should soon understand better the properties of the hot hadron 
fireball and its presumably deconfined structure. 
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